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Abstract 

The errors in the arguments, numerical results, and conclusions in the paper "Test of a modified 
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I. INTRODUCTION 



The modified BCS (MBCS) theory and its extension, the modified Hartree-Fock- 
Bogohubov (MHFB) theory, have been proposed and developed in a series of papers 
l,!^, y, y, 5| between 2001 and 2007. The aim of these theories is to restore the uni- 
tarity relation E? = R for the generalized single-particle density matrix i?, which is violated 
within the conventional BCS and HFB theories at finite temperature T. This restoration 
leads to an application of the secondary Bogoliubov transformation from the conventional 
quasiparticle operators to the modified quasiparticle ones. As a result the MBCS equation 
is obtained. The main merit of the MBCS theory is that it is a fully microscopic theory, 
which shows that the phase transition from the superfiuid phase to the normal one (the SN 
phase transition) is smoothed out in finite systems. The MBCS pairing gap does not collapse 
at the critical temperature T^. of the SN phase transition, where the conventional BCS gap 
vanishes. It also points out, for the first time, that the origin of the smoothing out of the 
SN-phase transition is the quasiparticle- number fiuctuations (QNF), which are ignored in 
the conventional BCS and HFB theories. The MBCS theory has been applied with success 
in several open-shell realistic nuclei and tested by its authors in an exactly solvable 

pairing model with Vt doubly folded levels % ^ , which is often referred to as the Richardson 
model, ladder model, or picket-fence model (PFM). Being an approximation that deals with 
the fiuctuations due to the system finiteness, the MBCS theory is sensitive to the size of the 
system under consideration. Within the half- filled PFM, e.g., i.e. when N = Vt, where is 
the number of particles, the MBCS gap decreases smoothly with increasing the temperature 
T up to a certain temperature Tm (> Tc), where a discontinuity occurs (See Refs. 
and references therein). For a single-particle spectrum with the level distance e equal to 1 
MeV, it has been found that his maximal temperature Tm decreases almost linearly with 
decreasing T from a value as high as Tm ~ 24 MeV for N = VL = 100 to a value as low as 
Tm — 0.7 MeV for N = Vt = Q In Ref. it has been pointed out, for the first time, 
that the asymmetry of the QNF with respect to the Fermi energy is the reason that causes 
the decrease of Tm with A^. In the same Ref. Q], it has also been found that it is sufficient 
to add just one more level, i.e. f2 = + 1, to restore the symmetry of the QNF up to much 
higher Tm- This allows to extend the temperature region of applicability of the MBCS 
theory to a temperature as high as T < 5 - 6 MeV for all values of the particle number 
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(for e = 1 MeV and the pairiiiEf interaction G = 0.4 MeV). The MBCS theory has been the 

nnn un 

subject of exchanging criticisms between us 15|, 16|] and the authors of Refs. P, 18|]. In a 



recent paper 



9|], the same authors reiterated their criticism by repeating the results of the 



same test within the PFM. Regrettably, beside largely duplicating the contents, which have 
already been published in Refs. , the authors of Ref. |9| added nothing 

new but further serious mistakes and false conclusions. The aim of the present article is to 

n 

point out the errors as well as the repetitions in Ref. [9(, which is referred to hereafter as 
the test y. 



II. INCORRECT STATEMENTS AND ERRONEOUS RESULTS 



A. Chemical potential 

The authors of the test [9] claim that the chemical potential of the PFM should not 
change with temperature T, remaining always in the middle of the single-particle spectrum 
taken to be zero in this case. This statement is true only within the BCS theory for the 
half-filled case (i.e. N = Q) if the self-energy correction term ~ Gvj is neglected in the 
single-particle energies. Within the BCS theory at T 7^ 0, the chemical potential A is 
determined as a Lagrangian multiplier so that the average particle number within the grand 
canonical ensemble (GCE) has the correct value. This is because the BCS theory violates the 
particle number, causing quantal fluctuations of particle number, and ignores the effects due 
to the QNF at finite temperature as well. In any other theory, including the MBCS theory, 
where these effects and/or correlations beyond the mean field are partially or fully taken into 
account, the chemical potential A becomes a function of T. An example is the well-known 
Lipkin-Nogami (LN) method 10|], which is an approximate particle- number projection to 
partially eliminate the particle-number fluctuations in the BCS theory. The LN method 
yields A = Ai + 2X2{N + 1). These chemical potentials are shown in Fig. [T]as functions of 
T for the PFM with N = Q = 10 and G = 0.4 MeV. So long as the pairing gap A exists, 
all of them change with T, and none of them remains at zero. 

The exact definition of the chemical potential fi is the difference between the ground-state 
(internal) energies of the systems with ± 2 particles HI], namely 



^[SiN + 2) - SiN)] , /i_ = ^[SiN) - 8{N - 2)] , 



(1) 
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FIG. 1: (Color online) Chemical potentials A, Ai and A2 as functions of T within the Lipkin-Nogami 
method for the PFM with N = Q = W and G = 0.4. 
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FIG. 2: (Color online) Chemical potentials as a functions of T at G = 0.4 MeV (a), and a 
function of G at T = (b) for the PFM with O = + 1 = 11 (a), and 17 = iV = 10 (b). The exact 
CE results are shown as the solid lines. In (a) the exact CE result for 17 = 10 is also shown as 
the thin solid line. The BCS and MBCS results are shown in (a) as the dotted and dashed lines, 
respectively. 

/i = i(/i+ + /i_) = ^[SiN + 2) - SiN - 2)] , (2) 

where S{N) denotes the ground-state energy of the A^-particle system in the zero- 
temperature case, or its internal (total) energy at finite temperature T. In the latter case, 
the thermal average within the GCE is carried out for a system which exchanges its energy 
and particle number with a heat bath. If the particle number of the system is fixed, the 
canonical ensemble (CE) should be used for thermal averaging. The PFM is a system with 
a fixed number of particles, where no particle number fluctuations are possible. Therefore 
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FIG. 3: (Color online) Entropies as functions of T for iV = O = 10 at G = 0.4 MeV. The 
thermodynamic entropies Sth are shown in panels (a), while panels (b) display the quasiparticle 
entropies Sqp. In the upper panels, (al) and (bl), these quantities are shown in the logarithmic 
scale, whereas the linear scale is used in the lower panels, (a2) and (b2). The notations for the 
lines are shown in (al), including the "exact results" of the test ^] (circles). In (bl) and (b2) 
the exact GCE and CE single-particle entropies are also shown as faint thick and thin solid lines, 
respectively. 

the meaningful exact results should be those obtained within the CE by using the partition 
function constructed from the exact eigenvalues of the pairing Hamiltonian [See Eqs. ([3]) for 
the expressions of the total energies within the GCE and CE]. As shown by the solid lines 
in Fig. O the chemical potential fi strongly depends on temperature T and the interaction 
parameter G. From Fig. [2] (a) one can also see that the performance of the MBCS theory 
(the dashed line) is quite reasonable. 
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B. Entropies 



Regarding the entropies, the results that the authors claim as the exact thermodynamic 
entropy, denoted as circles in Fig. 7 (a) of Ref. j^, are incorrect. The same circles are also 
shown in Fig. [3] (al). The correct exact GCE and CE results for this quantity are shown as 
thick and thin solid lines, respectively, in Fig. [3] (al) (in the logarithmic scale) and Fig. [3] 
(a2) (in the linear scale). They are obtained from the same Eq. (17) of the test [9,] with the 
total energies calculated within the GCE and CE by using the textbook definitions, namely 

with the grand partition function Z{(3, A), and partition function Z{(3, n) defined as 

Z(/3, A) = ^ e^^"Z(/5, n) , Z(/5, n) = ^p.(/5, n) , p,(/?, n) = d^^^e-^''"' . (4) 



In Eq. (jll), are the eigenvalues of the pairing Hamiltonian for the system with n particles, 
whereas (ii"'' are their degeneracies. The GCE sum in Eq. (JSj) is carried out over all n = 
1, . . ., Q-1 with blocking properly taken into account for odd n. As compared to these 
exact GCE values [thick solid line in Fig. [3] (al)], the results of the test ^ (circles) at T < 
0.1 MeV are wrong by more than one order of magnitude. Meanwhile, the thermodynamic 
entropy obtained within the MBCS theory (the dashed lines) are found sandwiched between 
the exact GCE and CE results. The quasiparticle entropy obtained within the MBCS theory 
agrees quite well with the BCS one, except for the region around the critical temperature Tc, 
where the BCS gap collapses. At high T they both converge to the exact GCE result for the 
single-particle entropy as expected [See faint thick solid hues Figs. [3] (bl) and [3] (b2)]. The 
exact GCE and CE single-particle entropies are obtained by using the occupation numbers 
fj on the jth single-particle orbital within the GCE and CE, respectively. The latter are 

(s) 

the ensemble averages of the exact state-dependent occupation numbers /j , namely 



with 



where (C^*"*)^ determine the weights of the eigenvector components, and Nj^"* are the partial 
occupation numbers weighted over the basis states |A;). 



III. REPETITIONS OF PUBLISHED MATERIAL 



Sections 2 and 3 of the test 



9| repeat the arguments of Ref. 



where the same authors 



rephed to the our comments in Ref. 



a]. This is clearly seen, e.g., in Fig. 1 of the test 



which is obtained by putting all the lines from the upper panels (for the pairing gaps) of 
Fig. 1 in Ref. Sj into one panel. The only difference is that Fig. 1 of Ref. 8|] is for N = 10, 
while the same figure in the test [s] is presented for N = 14. This modification obviously 
offers no new physics as compared to Fig. 1 in the previously published Ref. jsj. 

Section 3 of the test JOlalso repeats the discussion of Refs. . As we have 

pointed out in Refs. [3,0, ^, 6], the validity of the MBCS depends on how the QNF is taken 
into account. If the single-particle spectrum is too small, the QNF can be large even at the 
top and bottom of the spectrum, and its profile becomes strongly asymmetric with respect 
to the middle of the spectrum already a low T. This leads to some artifacts as negative 
values of the pairing gap, or a discontinuity in its temperature dependence. Therefore, 
a criterion is introduced to include the QNF symmetrically from both sides of the Fermi 
energy in such a way that the negative and positive wings discussed in Fig. 2 (b) of Ref. 

nearly cancel each other. For the PFM under consideration, it turns out that, for small 
A^, e.g. < 14 with the level distance equal to 1 MeV and G= 0.4 MeV, it is sufficient to 
enlarge the space by just one level = + 1 to achieve a symmetric profile of the QNF 
up to a rather high temperature, while adding or subtracting more levels again decreases 
the value of the limiting temperature, at which the QNF profile becomes asymmetric. For 
realistic heavy nuclei, such as ^^°Sn e.g., this criterion is often unnecessary so long as the 



In this respect. Figs. 2-4 



whole single^particle space is included (Fig. 1 - f of Ref. 
of the test 9(] and the related discussions therein are just a disguised tautology since they 
demonstrate nothing but the asymmetry of the QNF, which has already been shown and 



discussed previously in Fig. 2 of Ref. 



In Fig. 2 (b) of our paper j^, the partial components of the thermal gap are shown to 
clarify the temperature dependence of the total gap, which is a physical quantity. It seems 
that the authors of the test [q!] misunderstood this idea when they borrowed our results to 
divide the MBCS thermal gap into the hole and particle gaps [See Eqs. ^^^) (16) of 
the test [9i|]]- They summed up the values of 6Aj discussed in our paper [31] over the levels 
below the Fermi energy separately from those above it, and considered these two sums as 



^, and Fig. 1 of Ref. 
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the hole and particle gaps, respectively. Obviously, such separation of the thermal gap into 
the hole and particle gaps as two physical quantities is misleading. The MBCS thermal gap 
is level- independent, and induced by all the QNF, which should be included by summation 



over all quasiparticle orbitals. Figures 5 (c) and 5 (d) of the test jo] repeat the upper panel 
of Fig. 3 in Ref . [7| , which has previously been published by the same authors under almost 
the same title. Two additional panels (a) and (b) bear no additional physics information. 



These issues have been refuted by us in Ref. {g]. 

Last but not least, that the MBCS theory does not cure the particle-number fluctuation 
(PNF) has already been discussed long ago in Ref. where the PNF was calculated for 



several neutron-rich Ni isotopes [2]. This is the reason why several methods of particle- 
number projection were applied within the MBCS theory in Ref. |5|], whose calculations 
were carried out not only for the PFM but also for the realistic nucleus ^^°Sn. Therefore, 
Fig. 6 and the related discussions in the test [9|] are redundant because they just reiterate 
a well-known fact by using an oversimplified toy model with just two levels. 

IV. INTRODUCING MISLEADING QUANTITIES AND MAKING WRONG 
CONCLUSIONS 

A part from the errors and repetitions mentioned above, the authors of the test [9] also 
introduced misleading discussions by comparing the MBCS predictions of quasiparticle en- 
ergies with the "exact" ones (See Fig. 3 of the test [9]), whose definition remains completely 
obscure. A quasiparticle state is an approximation, which arises from the canonical Bogoli- 
ubov transformation, and is related to the mean-field concept. So long there is no exact 
mean field, there cannot be exact quasiparticle energies. The eigenvalues of the PFM, which 
are obtained by diagonalizing the pairing Hamiltonian or solving the Richardson's equa- 
tions, are not quasiparticle energies. In order to have a better agreement with the exact 
eigenvalues one needs to go beyond the conventional BCS theory to take into account the 
correlations beyond the mean field. It is for this purpose that the particle-particle self- 
consistent RPA (SCRPA) I12I, [3 1 (for low value of the pairing interaction parameter G) 
and quasiparticle RPA (SCQRPA) [iJ] (for an arbitrary value of G) were developed. But 
while the SCRPA and SCQRPA can describe rather well the low excitation energies, they 
cannot describe the complicate splitting of high-lying levels due to configuration mixing in 
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FIG. 4: (Color online) MBCS gaps A as functions of temperature T obtained within the PFM, 
which has = 10 particles, and Q = N + k levels, (a) Results obtained with G = 0.6 MeV, k = 
0, 1,. . . , 8; (b) Results obtained with G = 0.55 MeV, A; = 0, 1, 2, and 3. The numbers at the right 
margins of (a) and (b) show the corresponding values of O. The single-particle energies for the sets 
(a) and (b) are specified in the text. 



the exact solutions 



15|- 



Apart from the conclusions of the test |9'], which obviously become invalid as a conse- 
quence of the above analysis, the authors of the test [q^] also claim: " We confirm that there 
exists a single example of the PFM (with the number of levels Q equal to the number of 
particles N plus one ) in which the MBCS produces the thermal behavior of the pairing gap 
similar to the one of a macroscopic theory up to rather high temperatures. On the other 
hand, we demonstrate that in all other examples of the PFM with Vt ^ N the theory predicts 
phase transitions of unknown types at a much lower temperature^ Of course, this statement 
is incorrect. It is valid only for a PFM with = 10, G = 0.4 MeV with the level distance 
e = 1 MeV. It is not difficult to find counter examples, some of which are shown in Fig. HI 
It depicts the MBCS pairing gaps A obtained for the PFM with = 10, and VL = N + k, 
which are grouped in two sets. The set in Fig. H] (a) consists of nine MBCS gaps, which are 
obtained with A; = 0, . . . , 8 by using G = 0.6 MeV, and the following single-particle energies: 
€j = 2.1x j (MeV) for j < 3 and j > 7, whereas ej = es + j - 3 (MeV) for 4< j < 6. The 
set in Fig. H] (b) represents four MBCS gaps, which are obtained with A; = 0, 1, 2, and 3 by 
using G = 0.55 MeV, and the following single-particle energies: ej = 2.1 x (j — 5.5) (MeV) 
for j < 3 and 7 < J < 10, ej = j- 5.5 (MeV) for 4 < j < 6, and ej = 2.1 x (0.55 x eio +j - 10) 
(MeV) for 11 < j < 13. All thirteen MBCS gaps in Fig. Hlare smooth functions of T at 
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< T < 5 MeV. This figure clearly shows how changing the size of the system affects the tail 
of the gap at T > 2 MeV in these cases. These examples and the analysis in the present 
article are more than sufficient to rule out all the conclusions of the test . 

V. CONCLUSION 

In conclusion, the test not only repeats already published results, but also contains 
incorrect statements and detectable errors, which make its conclusions invalid. 

Acknowledgment s 

I thank N. Quang Hung for critically reading the manuscript and assistance in numerical 
calculations, which were carried out by using the Fortran IMSL Library by Visual Numerics 
on the RIKEN Super Combined Cluster (RSCC) system. 



[1] N. Dinh Dang and V. Zelevinsky, Phys. Rev. C 64, 064309 (2001); Ibid. 65, 069903(E) (2002). 

[2] N. Dinh Dang and A. Arima, Phys. Rev. C 67, 014304 (2003) 181; Ibid 68, 039902(E) (2003). 

[3] N. Dinh Dang and A. Arima. Phys. Rev. C 68, 014318 (2003). 

[4] N. Dinh Dang, Nucl. Phys. A 784, 147 (2007). 

[5] N. Dinh Dang, Phys. Rev. C 76, 064320 (2007). 

[6] N. Dinh Dang and A. Arima, Phys. Rev. C 74, 059801 (2006). 

[7] V.Yu. Ponomarev and A.I. Vdovin, Phys. Rev. C 72, 034309 (2005). 

[8] V.Yu. Ponomarev and A.I. Vdovin, Phys. Rev. C 74, 059802 (2006). 

[9] V.Yu. Ponomarev and A.I. Vdovin, Nucl. Phys. A 822, 1 (2009). 

[10] H.J. Lipkin, Ann. Phys. 9, 272 (1960); Y. Nogami, Phys. Rev. 134, B313 (1964); H.C. Prad- 

han, Y. Nogami, and J. Law, Nucl. Phys. A 201, 357 (1973). 

[11] P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer, Heidelberg, 2004). 

[12] J. Dukelsky and P. Schuck, Phys. Lett. B464. 164 (1999). 

[13] N.D. Dang, Phys. Rev. C 71, 024302 (2005); N.D. Dang, Phys. Rev. C 74, 024318 (2006); 
N.D. Dang and K. Tanabe, Phys. Rev. C 74, 024326 (2006). 



10 



[14] N.Q. Hung and N.D. Dang, Phys. Rev. C 76, 054302 (2007); Ibid. 77, 029905(E) (2008); N.D. 

Dang and N.Q. Hung, Phys. Rev. C 77, 064315 (2008). 
[15] E.A. Yuzbashyan, A. A. Baytin, and B.L. Altshuler, Phys. Rev. B 68, 214509 (2003). 



11 



